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Abstract 

When nonlinear effects on tfie gluon evolution are included with constraints from 
HERA, the gluon distribution in the free proton is enhanced at low momentum frac- 
tions, X ^ 0.01, and low scales, ^ 10 GeV^, relative to standard, DGLAP-evolved, 
gluon distributions. Consequently, such gluon distributions can enhance charm pro- 
duction in pp collisions at center of mass energy 14 TeV by up to a factor of five at 
midrapidity, y ~ 0, and transverse momentum p-r in the most optimistic case. We 
show that most of this enhancement survives hadronization into D mesons. Assuming 
the same enhancement at leading and next-to-leading order, we show that the D en- 
hancement may be measured by reconstruction in the K'lr'^ decay channel with 
the ALICE detector. 
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1 Introduction 



The parton distribution functions, PDFs, of the free proton are determined through 
global fits obtained using the leading-order, LO, next-to-leading order, NLO, or even 
next-to-next-to-leading order, NNLO, formulation of the Dokshitzer, Gribov, Lipatov, 
Altarelli and Parisi, DGLAP, scale evolution equations [Tj. In particular, the HERA 
data on the proton structure function F2{x,Q'^) |2j as a function of Bjorken-x and 
squared momentum transfer Q^, and, especially, the slope, dF2{x, Q"^) /dlnQ'^ , in 
the small-x, 3 x 10"^ < x < 5 x 10"^, and small-Q^ region, 1.5 < ~ 10 GeV^, set rather 
stringent constraints on the small-x gluon distributions. The agreement of the global 
fits with the measured F2{x,Q'^) is, in general, very good but certain problems arise. 
When the small-x and small-Q^ region is included in the DGLAP fits, they are not as 
good as the excellent ones obtained at larger values of x and [H]- In addition, some 
NLO gluon distributions |H become negative at small x for on the order of a few 
GeV\ 

The kernels of the DGLAP equations only describe splitting of one parton into two 
or more so that the resulting equations are linear in the PDFs. This ignores the fact 
that, at low Q^, the small-x gluon density may increase to the point where gluon fusion 
becomes significant. These fusions generate nonlinearities in the evolution equations. 
The first nonlinear corrections, the GLRMQ terms, were derived by Gribov, Levin 
and Ryskin and also by Mueller and Qiu [H]. Eventually, at even smaller x and Q^, 
nonlinearities are expected to dominate the evolution to all orders. This fully nonlinear 
region, where both the linear DGLAP evolution and the GLRMQ-corrected DGLAP 
evolution are inapplicable, is the gluon saturation region, see e.g. Ref. |n|. 

Outside the saturation region, incorporating the nonlinearities may improve the 
global fits when the small- a; and regions are included. Recent work in Ref. jlj, where 
the LO DGLAP evolution equations were supplemented by the GLRMQ terms, showed 
that the nonlinearly-evolved PDFs reproduce the HERA F2 measurements at x > 3 x 
10"^ and Q2 > 1.5 GeV^ equally well or even better than the conventional LO PDFs 
such as CTEQ6L 0. The nonlinearly-evolved gluon distributions at ^ 10 GeV^ and 
X ^ 0.01, however, were clearly enhanced relative to CTEQ6L and CTEQ61L P|. As 
shown in Fig. 1 of Ref. jini, the enhancement arises because the nonlinear evolution is 
slower than DGLAP alone. At higher x and the nonlinear and linear evolution of the 
gluon distributions should become very similar to fit the same data. An enhancement 
can also be expected at NLO. However, since the NLO small-x gluon distributions are 
typically reduced relative to LO, at NLO the enhancement may be smaller than at 
LO [3]. 

Since the same HERA data can be reproduced by linear evolution starting from 
a relatively fiat gluon distribution and by nonlinear evolution with clearly enhanced 
small-x gluons, other observables are necessary to probe the effects of the nonlinearities. 
In Ref. charm production in pp collisions at the LHC was suggested as a promising 
candidate process. Due to gluon dominance of charm production and the small values 
of X and Q"^ probed, x ^ 2 x 10"*^ and Q'^ ^ 1.69 — 6 GeV^ at midrapidity and 
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transverse momentum^ ~ 0, charm production at the LHC is sensitive to the gluon 
enhancement. The resulting charm enhancement was quantified in Ref. [in| by the LO 
ratios of the diff'erential cross sections computed with the nonlinearly-evolved EHKQS 
PDFs [H, obtained from DGLAP+GLRMQ evolution, relative to the DGLAP-evolved 
CTEQ61L PDFs. 

The enhancement of the nonlinearly-evolved gluons increases as x and Q"^ decrease. 
Consequently, the charm enhancement increases with center of mass energy, y/s. Thus 
the maximum enhancement at the LHC will be at a/s = 14 TeV and small charm 
quark transverse momentum. The sensitivity of the charm enhancement to the value 
of the charm quark mass, mc, as well as to the choice of the factorization, and 
renormalization, Q^, scales was studied in Ref. assuming — Q\ — Q\, ^ Tfiq^, 
the charm transverse mass squared, mf. = + ml. For the most significant charm 
enhancement, rUc and Q'^/m'^ should both be small. A comparison of the NLO to- 
tal cross sections with low energy data shows that the data prefer such small rric and 

combinations [IHII21- The smallest scales and thus the largest enhancement are 
obtained with rric = 1.3 GeV and = m\. In this case, the ratio of the inclusive differ- 
ential cross section, (Pa/dpTdydy2, computed with EHKQS set 1 relative to CTEQ61L 
is greater than 5 for rapidities \y,y2\ ~ 2 where y and y2 are the c and c rapidities, 
respectively. 

In Ref. [E], the enhancement was described only for charm production. Neither its 
subsequent hadronization to D mesons nor its decay and detection were considered. 
In this paper, we address these issues to determine whether the charm enhancement 
survives hadronization and D decay. At the LHC, the ALICE detector |Tni is per- 
haps in the best position for measuring such an enhancement since it is capable of 
reconstructing hadronic decays down to very low transverse momentum. 

We first consider how much of the LO charm enhancement survives in the final 
state D meson distributions. Charm quarks are hadronized using the PYTHIA string 
fragmentation model |14|. We show that, for the most optimistic case with a factor of 
five charm enhancement for pt ^ 0, the D enhancement is a factor of three for — 0. 

Since the ALICE detector allows direct measurement of the D meson pt distribu- 
tion through reconstruction in the K~7r~^ decay channel, we then determine whether 
or not the surviving D enhancement can be detected above the expected experimental 
statistical and systematic uncertainties. To determine realistic statistical uncertainties, 
we calculate the NLO cross section in the way most compatible with our LO enhance- 
ment, as described below. Then, using the error analysis developed by one of us (A.D.) 
in Ref. [IH], we demonstrate that detection of the enhancement is possible. 

Finally, we consider whether NLO charm cross sections, calculated with linearly- 
evolved PDFs and different combinations of mc, Qp and Q/j, can mimic the charm 
enhancement. Our results show that this is unlikely. 

This paper is organized as follows. In Section 2, we describe our charm calculations 
and define how the NLO cross section most compatible with the LO enhancement is 

^Here we use pr for the transverse momentum of the charm quark and for the transverse 
momentum of the D meson. 
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computed. Hadronization and reconstruction of D mesons are considered in Sections 
3 and 4, respectively, along with a discussion of the experimental uncertainties. We 
then generate "data" based on the enhanced cross sections and the experimental un- 
certainties. These data are then compared to compatible NLO calculations to learn 
whether the enhancement is measurable for a unique set of parameters in Section 5. 
We conclude in Section 6. 

2 Charm enhancement from nonhnear PDF evolution 

According to collinear factorization, the inclusive differential charm hadroproduction 
cross sections at high energies can be written as 

dapp^ccxiVs,mc,Q'^ji,Qp) = J2 M^i^Qf) fji^2,Q%) 

i,j=q,q,9 

(g) dcrij^cc{k}{as{Q%), Q%, rUc, Xi, X2), (1) 

where daij^cc{k} is the perturbative partonic hard part, calculable as a power series 
in the strong coupling agiQ^). The proton PDFs for each parton at fractional 
momentum Xi(x2) and factorization scale Qp are denoted by fi{x,Qp). At LO, where 
da oc a;^(Q|.), only the subprocesses gg cc and qq cc are allowed [H] so that 
{k} = 0. At NLO, where da oc subprocesses where {k} 7^ 0, e.g. gg — > ccg 

and gq — > ccq contribute. The gq channel, new at NLO, only contributes a few percent 
of the total cross section. 

The charm production enhancement studied here and in Ref. |lfl| results from 
the nonlinearly-evolved EHKQS PDFs^ where the gluon distribution is enhanced for 
X ^ 0.01 at the few-GeV scales. The EHKQS PDFs were constructed in Ref. [I| us- 
ing CTEQ5L [m and CTEQ6L as baselines with the HERA data as constraints. 
The EHKQS sets have initial scale Ql = 1.4 GeV^ and a four- flavor Aqcd value of 
^QCD — 0.192 GeV. Following Ref. [101, we quantify the charm enhancement against 
charm production computed with the CTEQ61L LO PDFs where the data were fit with 
the one-loop a,. The CTEQ61L set takes Ql = 1.69 GeV^ and AgJ^^ = 0-215 GeV. For 
consistency, we calculate at one loop with the appropriate value of Aq^^q for each 
set. 

Previously [10], we worked at LO only since the EHKQS sets are evolved according 
to the LO DGLAP+GLRMQ equations using a one-loop evaluation of a^. Thus these 
LO distributions should generally not be mixed with NLO matrix elements and the 
two-loop as- However, the charm quark total cross section is increased and the pt 
distribution is broadened at NLO relative to LO [iH]. Thus, to determine whether 
or not the enhancement is experimentally measurable, we must go beyond the ratio 
presented in Ref. [10]. To accomplish this, we assume that the enhancement will be 
the same at NLO as at LO and employ a NLO cross section closest to the calculation 
of the enhancement in Ref. P)J. 

^These PDFs are available at www.urhic.phys.jyu.fi. 
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As described in Ref. [TH], the theoretical K factor may be defined in more than one 
way, depending on how the LO contribution to the cross section is calculated. In all 
cases, the C(a^) contribution to cross section is calculated using NLO PDFs and the 
two-loop evaluation of a^. If the LO contribution is also calculated using NLO PDFs 
and a two-loop a^, this is the "standard NLO" cross section. It is used in most NLO 
codes, both in the global analyses of the NLO PDFs and in evaluations of cross sections 
and rates |TH]. The K factor formed when taking the ratio of the "standard NLO" cross 
section to the LO cross section with the NLO PDFs [IS.], -f^o^\ indicates the convergence 
of terms in a fixed-order calculation [121. On the other hand, if the LO contribution 
to the total NLO cross section employs LO PDFs and the one-loop a^, we have a cross 
section which we refer to here as "alternative NLO". The K factor calculated taking 
the ratio of the "alternative NLO" cross section to the LO cross section with LO PDFs 
[T8] . \ indicates the convergence of the hadronic cross section toward a result. If 
K^f^ > K2'\ convergence of the hadronic cross section is more likely ^O]- This is indeed 
the case for charm production jTH|. We also note that K2'^ is a much weaker function 
of energy than Kq^\ Since, in the absence of nonlinear NLO PDFs, the "alternative 
NLO" cross section is more consistent with the enhancement calculated in Ref. [Tfl] . 
we use this cross section to calculate the NLO D meson rates and pt spectra. We note 
also that, in both cases, the pt distributions have the same slope even though K2'\ for 
the alternative NLO cross section, is somewhat smaller. Thus, using a non-standard 
NLO calculation will not change the slope of the pr distributions, distorting the result. 

The LO and NLO calculations used to obtain the full NLO result in both cases can 
be defined by modification of Eq. Jl]). For simplicity, we drop the dependence of the 
cross section on a/s, rric, Q% and Q|. on the left-hand side of Eq. (|T} in the following. 
We thus define the full LO charm production cross section as 

where the superscript "LO" on daij^cc indicates the use of the LO matrix elements 
while the superscript "IL" indicates that the one-loop expression of as is used. The LO 
cross section typically used in NLO codes employs the NLO PDFs and the two-loop 
(2L) as so that 

In either case, the NLO contribution, 0{a^g) for heavy quark production, is 

where the superscript "NLO" on ddij^cck indicates the use of the NLO matrix elements. 
The additional sum over k in Eq. ^ includes the virtual, /c = 0, and real, k = q, q 
or g depending on i and j, NLO corrections. In the calculations of da^ and dao{a3), 
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we use the value of Aq^j-, given for the NLO PDFs and work in the MS scheme. The 
standard NLO cross section is then 

'^'^NLO = dcr^ + dao{ai) (5) 
while our "alternative NLO" cross section is defined as 

do'NLO = d-'^LO + d(Jo(^a3) ■ (6) 

Since the enhancement in Ref. [E] was defined using daj^ only, the best we can do is 
to use the alternative NLO cross section in our analysis, as described below. 

We now discuss how the enhancement is taken into account in the context of the 
NLO computation. We calculate the LO inclusive charm pt distribution, d'^a/dpj'dy, 
with the detected charm (anticharm) quark in the rapidity interval Ay with \y\ < 1, 
motivated by the pseudorapidity acceptance of the ALICE tracking barrel, \ri\ < 0.9. 
The rapidity, j/2, of the undetected anticharm (charm) quark is integrated over. The 
charm enhancement factor R{pt, Ay) is then 



dy J dy2^ 
dy j dy2 



Cr^EHKQSj 

Ay"'''' J """" dpTdydy2 

(t(cTEQ61LI 



RiPT. Ay) = : i^LlZU ■ (7) 

lAy""" J ""^ dpTdydy2 



Numerically, this ratio is very close to R{pT,y,y2), computed in Ref. as seen by a 
comparison of R{pT,Ay) in Fig. [T] with Fig. 2 of Ref. [TO] . 

Next, we assume that the enhancement calculated at LO is the same when calcu- 
lated at NLO. This is a rather strong assumption but, until the nonlinear evolution 
has been completely analyzed to NLO, it is the only reasonable assumption we can 
make to test whether the enhancement can be detected with ALICE which will mea- 
sure the physical distribution. The alternative NLO cross section is therefore the 
closest in spirit to the LO computation in Ref. [IH]. Thus, the enhanced NLO charm 
Pt distribution is 

R{pT,Ay)dcrt^oiAy)/dpT. (8) 

In our calculations, we use values of the charm quark mass and scale that have been 
fit to the total cross section data using standard NLO calculations. The best agreement 
with the total cross section data is obtained with rric = 1.2 GeV and = 4m^ for 
DGLAP-evolved NLO PDFs such as CTEQ6M and MRST Nearly equivalent 
agreement may be obtained with nic = 1.3 GeV and = ml [II1II2]. Agreement with 
the fixed-target total cross sections can only be achieved with higher rric by making 
the factorization scale, Qp, larger than the renormalization scale, Using a lower 
value of increases the cross section by infiating a^. If Qp < Ql, the PDFs are 
unconstrained in and are thus unreliable. We keep Qp = Q\ since all typical PDFs 
are fit using this assumption. Thus we limit ourselves to relatively small values of rric 
to obtain agreement with the total cross section data. 
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We note that while rric is the only relevant scale in the total cross section, is 
used instead of in the calculations of and da^l^Q^Ay) / dpr to control py-dependent 
logarithms at NLO [101. Our main results are then based on the inputs that give the 
best agreement with the total cross section data, rric = 1.2 GeV and = Am\ as 
well as mc = 1.3 GeV and = m^. These two choices will form the baseline results 
against which other parameter choices will be compared to see if the enhancement can 
be detected. 



3 Prom charm to D enhancement 



Previously ^Hj, we did not include parton intrinsic transverse momentum, fcy, broaden- 
ing or fragmentation. Since the effect of intrinsic kx is quite small at LHC energies, on 
the order of 10% or less [TT], we have not included intrinsic kx in our calculations. To 
make a more realistic D meson distribution, we have modified the charm pt distribu- 
tion by the heavy quark string fragmentation in PYTHIA [T^, as explained below. The 
resulting D distribution is significantly harder than that obtained using the Peterson 
fragmentation function |21| . 

We first show how the pT-dependent enhancement, calculated for the charm quark, 
is refiected in the D meson pr distribution. Charm events in pp collisions at ^/s = 
14 TeV are generated using PYTHIA (default settings) with the requirement that one 
of the quarks is in the interval \y\ < 1. The charm quarks are hadronized using the 
default string model. Since c and c quarks fragment to D and D mesons'^, respectively, 
in each event related {c,D) and (c, D) pairs can easily be identified^. These pairs 
are reweighted to match an arbitrary NLO charm quark pt distribution, dN^-^Q/dpr- 
If dNpYTniAl dpT is the charm pt distribution given by PYTHIA, each (c, D) pair is 
assigned the weight 

where Pt is the transverse momentum of the charm quark of the pair. Therefore, the 
reweighted final-state D distribution corresponds to the one that would be obtained 
by applying string fragmentation to the NLO c-quark distribution. 

In Fig. ^ we compare the enhancement factor i?, calculated in Eq. for c quarks 
and D mesons generated from the weighted PYTHIA charm distributions. The two 
cases described previously, rric = 1-2 GeV, = Arri^ (left-hand side) and rric = 
1.3 GeV, = rri^ (right-hand side) are considered. In both cases, the enhancement 
survives after fragmentation. It is interesting to note that the D enhancement is 
somewhat lower than that of the charm: in the most optimistic case, the factor of five 
charm enhancement has reduced to a factor of three for the D mesons. This occurs 
because, for a given p^, the D spectrum receives contributions from charm quarks with 

•^Here D = D+,D°. 

"^Events containing charm baryons were rejected. 
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Figure 1: Enhancement factor R{pT-,^v) for charm quarks (dashed histogram) and for 
D (= D~^,D^) mesons (sohd histogram), obtained after PYTHIA string fragmentation. The 
left-hand side shows the result for rric = 1.2 GeV and = Arrij^ while the right-hand side is 
the result for rric = 1-3 GeV and = m^. 



Pt ~ Pti where the charm enhancement is smaller. The D enhancement also vanishes 
with increasing transverse momenta, like the charm enhancement. 



4 reconstruction in pp collisions with ALICE 

The transverse momentum distribution of mesons produced at central rapidity, \y\ < 
1, can be directly measured from the exclusive reconstruction of D° K^n^ decays 
(and charge conjugates) in the Inner Tracking System (ITS), Time Projection Chamber 
(TPC) and Time Of Flight (TOF) detectors of the ALICE barrel, |r/| < 0.9 [IS|. The 
main feature of the -D° decay topology is the presence of two tracks displaced from the 
interaction point by, on average, 50 /xm, for ~ 0.5 GeV, to 120 /xm, for 5 GeV. 
Such displacement can be resolved with the ALICE tracking detectors and thus a large 
fraction of the combinatorial background in the K^tt"^ invariant mass distribution can 
be rejected. The low value of the magnetic field, 0.4 T, and the K/n separation in 
the TOF detector extend the D° measurement down to ~ 0. The analysis strategy 
and the pertinent selection cuts were studied with a realistic, detailed simulation of 
the detector geometry and response, including the main background sources [T5l l22] . 

The expected ALICE performance for pp collisions at ^/s = 14 TeV is summarized 
in Fig. El where the estimated relative uncertainties are reported as a function of p^. 
The main contributions to the py-dependent systematic error (triangles) are the de- 
tector acceptance and reconstruction efficiency corrections (squares), ~ 10%, and the 
correction for feed-down from bottom decays, B ^ + X (open circles), ~ 8%. The 
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Figure 2: Estimated relative uncertainties on the measurement of the differential cross 
section in pp collisions at the LHC with ALICE jTH]. Statistical uncertainties correspond to 
10^ minimum-bias pp events (an ~ 9 month run with a luminosity of ~ 5 x 10^'' cm^^s^^). 



latter is estimated based on the present 70 — 80% theoretical uncertainty in the bb cross 
section at LHC energies [231 ■ However, v^^e expect this uncertainty to be significantly 
reduced by the measurement of B decays to single electrons, i? — > + X, in AL- 
ICE [2HI. The pT-independent systematic error introduced by normalization to the pp 
inelastic cross section (inverted triangles) is also reported. This cross section v^^ill be 
measured by the TOTEM experiment [25 with an ~ 5% uncertainty. 

The statistical error corresponding to 10^ minimum-bias pp events (filled circles), 
an ~ 9 month run with a luminosity of ~ 5 x 10'^° cm~^s~^, is smaller than or on the 
order of the p^-dependent systematic error up to pf^ ~ 24 GeV for the alternative NLO 
cross section calculated using nic = 1.2 GeV, = Am^ and the CTEQ6 PDFs with 
no enhancement. The relative statistical error depends on the charm cross section, 
as we now explain. For a given p^ or p^ range, the statistical error is the error 
on the number of real (^) mesons in the K^tt^ invariant mass distribution, the 
signal, S(p^). The error is equal to ^S^pfp) + B(p^)/S(p^) where B(py) is the number 
of background candidates in the Z)° mass region. Then, at low p^, the error is ~ 
y^B(p^)/S(py) oc l/{dao/dp^) since the invariant mass distribution is dominated by 
combinatorial background. At high p^, the background is negligible and the error 
becomes ~ 1/ \J'^{Pt) ^/ \J da^l dp^^. In our subsequent results, the statistical errors 
are calculated taking this cross section dependence into account. 
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5 Sensitivity to the enhancement 



Figure ini shows the double-differential cross section, a o / dp^dy , in \y\ < 1 as 
a function of the transverse momentum. The points represent the expected "data" 
measured by ALICE, obtained from the alternative NLO cross section scaled by the 
enhancement factor R{pT,Ay) defined in Eq. and modified by string fragmen- 
tation. The solid and dashed curves are obtained by applying string fragmentation 
to the alternative NLO and standard NLO cc cross sections, respectively. Thus, the 
"data" points include the enhancement while the curves do not. The horizontal error 
bars indicate the bin width, the vertical error bars represent the statistical error and 
the shaded band gives the py-dependent systematic error. The 5% pT-independent 
systematic error on the normalization is not shown. The left-hand side shows the 
results for rric = 1.2 GeV and = 4m^ while the right-hand side shows those for 
rric = 1.3 GeV and = m\. The standard NLO cross section, Eq. and the Oi^a^g) 
contribution to the alternative NLO cross section, Eq. ((3)), were calculated using the 
HVQMNR code [2S1 with CTEQ6M and AgJ^o = 0-326 GeV. The LO contribution to 
the alternative NLO cross section, Eq. was calculated using the CTEQ61L PDFs. 
Fragmentation was included as described in Section 3. The enhancement, the differ- 
ence between the data and the solid curve visible for ^ 3 GeV, is more pronounced 
for the larger mass and lower scale, shown on the right-hand side of Fig. [HI 

There is a significant difference between the alternative and standard NLO distri- 
butions. Part of the difference is due to the one- and two-loop evaluations of as since 
of" < al^. This decrease will in turn reduce the 0{al) contribution to the alterna- 
tive NLO result relative to the LO component of Eq. (jS)). In addition, the standard 
NLO cross section would be reduced overall relative to a calculation with the same 
Aqcj-, at LO and NLO. However, these factors alone cannot explain the rather large 
difference between the standard and alternative NLO cross sections at low p^. The 
most important contribution is the large differences between the LO and NLO gluon 
distributions, especially at low scales. The slope of the CTEQ61L gluon distribution 
at = 1.69 GeV^ with x is very small until x > 0.01. On the other hand, the 
CTEQ6M gluon x slope is large and has the opposite sign relative to CTEQ61L for 
X < 0.04. The ratio of the two sets <it x ^ 10"^ is very large, CTEQ61L/CTEQ6M 
~ 100. At = 5.76 GeV^, the scale corresponding to 4m^ with nic = 1.2 GeV, this 
ratio decreases to a factor of two. We note that at fixed-target energies, ^/s < 40 GeV, 
the standard and alternative NLO results are indistinguishable from each other since 
the LO and NLO gluon distributions are rather similar in this relatively high x region, 
0.05 <x< 0.1. 

In order to address the question of the experimental sensitivity to the effect of 
nonlinear gluon evolution on low-pT charm production, we consider, as a function of 
Py, the ratio of the simulated data, including the enhancement, to alternative NLO 
calculations using a range of mc and along with PYTHIA string fragmentation. We 
denote this ratio as "Data/Theory". Thus, given the measured px distribution, we 
try to reproduce this result with NLO calculations employing recent linearly-evolved 
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Figure 3: Comparison of tlie simulated ALICE data generated from R{pt: ^y)d(y^i^o with tlie 
alternative (solid) and standard (dashed) NLO calculations. The effect of string fragmentation 
is included in the "data" points as well as in the curves. The left-hand side shows the result 
for rric = 1.2 GeV and = 4m^ while the right-hand side is the result for rric = 1.3 GeV 
and = rrij,. The error bars on the data represent the statistical error and the shaded band 
represents the pr-dependent systematic error. The 5% normalization error is not shown. 

PDFs and tuning rric and . We note that these parameters are not really free but are 
bounded by the range 1.2 < rric ^ 1-8 GeV and 1 ^ Q"^ /m\ < 4, as described in Section 
2 and in Ref. [El. 

Since the enhancement has disappeared for ^ 5 GeV, we refer to this unenhanced 
region as high p^. The region below 5 GeV, where the enhancement is important, 
is referred to as low p^. If no set of parameters can describe both the high- and low-p^ 
components of the distribution equally well, and, if the set that best reproduces the 
high-p^ part underestimates the low-p^ part, this would be a strong indication of the 
presence of nonlinear effects. 

The Data/Theory plots are shown in Fig.lH The points with the statistical (vertical 
bars) and p^-dependent systematic (shaded region) error correspond to the data of 
Fig. El including the enhancement, divided by themselves, depicting the sensitivity 
to the theory calculations. The black squares on the right-hand sides of the lines 
Data/Theory = 1 represent the 5% pT-independent error on the ratio coming from the 
cross section normalization. As clearly shown in Fig.O this error is, however, negligible 
with respect to the present estimates of the other systematic uncertainties (~ 13%). 
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Figure 4: Ratio of tlie generated ALICE data relative to calculations of the alternative NLO 
cross sections with several sets of parameters and PYTHIA string fragmentation. The left- 
hand side shows the result for uic = 1.2 GeV and = 4m while the right-hand side is the 
result for rric = 1.3 GeV and = m%. 



On the left-hand side, the thick solid curve with rric = 1-2 GeV and = 4my best 
agrees with the high-p^ ratio by construction since i? ~ 1 at large p^. It also shows the 
effect of the enhancement well beyond the error band for p;^ < 2 GeV. Better agreement 
with the data over the entire range can be achieved only by choosing a charm quark 
mass lower than 1.2 GeV, below the nominal range of charm masses, as illustrated by 
the dashed curve for mc = 1.1 GeV. Higher masses with = 4m\, produce much 
larger Data/Theory ratios than the input distribution. Choosing e.g. rric = 1-8 GeV 
(not shown) would give a larger Data/Theory ratio than the rric = 1.5 GeV result (dot- 
dashed curve). The ratio with rric = 1.3 GeV and = rri^ (dot-dot-dashed curve) 
gives a much larger ratio at low and drops below the data for > 8 GeV. 

We have checked how the results change when the renormalization and factorization 
scales are separated. When rric = 1-3 GeV, Qj^ = and Qp = Arri^, the faster 
evolution of the higher Qp and the larger as{Q%) resulting from the lower leads 
to reasonable agreement between data and theory at low p^. However, at high p^, the 
theory distribution is harder so that the Data/Theory ratio drops below the error band 
for pp > 2 GeV. On the other hand, when rric = 1.3 GeV, Q|. = 4m^ and Qp = m'^, 
the theory cross section is reduced relative to the data and the Data/Theory ratio is 
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above the error band over all p^. 

We also present the ratio using the MRST parton densities (MRST2001 LO H] in 
Eq. and MRST2002 NLO ^ in Eq. ^) with = 1.2 GeV and = 4^2^. vVe 
find that this result, the thin solid curve, also agrees reasonably well with the CTEQ6 
results shown in the thick solid curve for the same rric and . Thus, the enhancement 
seems to be rather independent of the PDF. The CTEQ61L and the MRST2001 LO 
distributions are similar at low x, suggesting that PDFs based on this MRST set would 
produce an enhancement like that of Ref. |^. However, the MRST2002 NLO and 
CTEQ6M NLO gluon distributions are very different at low x. The MRST2002 NLO 
gluon distribution is negative at low scales while the CTEQ6M gluon distribution goes 
to zero as X — *■ 0. Thus the effects of nonlinear evolution at NLO could be considerably 
different. 

On the right-hand side of Fig. |31 with = 1.3 GeV and = m^, the thick 
solid curve, employing the same parameters as the data, gives the best agreement at 
high p^. We note that even though the results with = Am^ and rric < 1.3 GeV lie 
closer to the data at low and within the combined statistical and systematic error at 
higher p^, the curves with these parameters have the wrong slopes for ^ 8 GeV. The 
systematic errors are more likely to shift the data points up or down as a whole rather 
than twist the shape. The statistical sensitivity is expected to be good enough 
to distinguish the difference in curvature. Varying Q']p and Q\ separately results in 
similarly poor agreement as that noted for mc = 1.2 GeV and = 4m^. Finally, the 
results obtained with the MRST PDFs, shown in the thin solid line, do not alter the 
conclusions already drawn for CTEQ6. 

6 Conclusions 

With constraints from HERA, the nonlinear DGLAP+GLRMQ evolution at LO leads 
to an enhancement of the free proton gluon distributions at a; < 0.01 and < 10 GeV^ 
relative to DGLAP-evolved LO sets such as CTEQ61L. Consequently, charm hadropro- 
duction at y/s ^ 1 TeV should be larger than expected from DGLAP-evolved PDFs 
alone [10]. In this paper, we have studied whether the EHKQS gluon distributions [Z| 
could generate an observable D meson enhancement in pp collisions at the LHC. Since 
larger x values are probed at lower energy colliders, the enhancement described here 
would be reduced. At RHIC, y/s = 200 GeV, the effect is too small to be reliably 
observed. However, D mesurements at the Tevatron, y/s = 1.96 TeV, may allow to 
detect an enhancement if the minimum p^ was lowerd to ~ 1 GeV. 

In order to consider more realistic p^ distributions and yields, we have calculated 
the NLO contribution to charm production using the HVQMNR code jSH]. Since the 
LO EHKQS PDFs cannot be used consistently with the NLO matrix elements, we 
assume the charm enhancement is the same at LO and NLO. We note that nonlinear 
effects on the NLO gluon distributions may be smaller than at LO, thus reducing the 
NLO charm enhancement. Therefore, our results may be considered upper limits of 
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the NLO D enhancement. Note also that if NLO DGLAP+GLRMQ PDFs that fit 
the small-x and small-Q^ HERA data were available, it would be possible to base our 
analysis on the standard NLO charm cross section instead of the "alternative NLO" 
result defined in Eq. Improved gluon distributions at low x and may make 
the standard and alternative NLO results more similar at high energies, as they are at 
lower ^/s where x is larger. 

Using the EHKQS LO PDFs and LO matrix elements for charm quark produc- 
tion and PYTHIA string fragmentation for D meson hadronization, we have demon- 
strated that more than half of the charm enhancement relative to calculations with the 
CTEQ61L LO PDFs indeed survives to the D mesons. In the most optimistic case, 
rric = 1.3 GeV and = m\, the factor of five charm enhancement at \y\ < 1 and 
Pt — ^ is reduced to a factor of three at 0. For larger values of and Q^, 

the charm enhancement is smaller because the gluon enhancement due to nonlinear 
evolution decreases with increasing . 

The D meson enhancement, however, drops rapidly with transverse momentum so 
that for p^^ 5 GeV it is only a few percent. Therefore, D measurement capability 
at small is necessary to verify the effect experimentally. The ALICE detector 
can do this through direct reconstruction in the K~7i~^ decay channel. We have 
demonstrated, using the error analysis of Ref. |lSj, that, in the most optimistic case, the 
enhancement can be detected above the experimental statistical and systematic errors. 
The sensitivity of the D enhancement to the scale has also been considered and we have 
shown that when the charm mass is somewhat smaller, rUc = 1.2 GeV, but the scale is 
larger, = Am^, it is more difficult to detect the enhancement over the experimental 
uncertainties. The ALICE sensitivity to D meson production at very low transverse 
momentum may further improve by combining the K^tt^ measurement with 

those of K^n^n^ and — > K~ti^ . A fast-simulation feasibility study of 

K^-n^-n^ reconstruction [23 indicates that a performance similar to that of 

— > i^r~7r+ could be achieved. More detailed analyses, currently in progress, will 
assess the low-p^ reach of this channel. 
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